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ABSTRACT 


The  present  extensive  use  of  the  ejector  as  a  fluid 
pump  has  brought  into  focus  the  need  for  an  adequate  ejector 
design  procedure.  An  analysis  is  presented  that  provides  a 
relation  for  ejector  behavior  in  terms  of  the  governing  non- 
dimensional  parameters.  Model  tests  are  presented  which 
a)  verify  the  analysis,  and  b)  establish  the  numerical  values 
for’  the  two  design  coefficients  that  are  required. 

The  intended  application  of  this  procedure  was  origin¬ 
ally  one  involving  the  exhaust  of  a  gas  turbine  as  the  driving 
fluid.  The  numerical  values  of  the  parameters  have  been 
confined  to  a  range  consistent  with  this  application,  and  in 
this  sense  only  the  analysis  is  limited. 

This  study  is  a  direct  extension  of  the  program  pre¬ 
viously  reported  in  references  (l)  and  (2). 


ACKNOWLEDGEMENT 

The  advice  and  suggestions  of  Professor  A.  L.  London 
are  most  sincerely  appreciated. 


iii 


TABLE  OF  CONTENTS 


rage 

vi 


References  . 

Nomenclature  . . vli 

Introduction  . .  1 

Summary  and  Conclusions . 4 

Analysis  ......  .  ,  ,  j 

Model  Test  Results  . . 19 

Application  to  Design . 29 

Appendix  A:  Velocity  and  Temperature  Traverse 

Results . . 35 

Appendix  B:  Analysis  of  Experimental 

Uncertainties . 4l 

Append!..  C:  Analysis  of  Pitot  Tube  Errors . 44 


LIST  OF  FIGURES 

Figure 

1.  Simple  Ejector  Without  a  Diffuser  .  18 

2.  Alternative  Locations  for  the  Primary  Jet  .  18 

3.  Ejector  Performance  -  Influence  of 
Departures  from  Simple  One-Dimensiona] 


Theory . 18 

4,5-  Ejector  Performance  -  Comparison  of 

Predicted  and  Model  Test  Results . 18 


6.  X  as  a  Function  of  W*  and  W*T*0*^.  .  .  25 

7.  Km  as  a  Function  of  W*T*° *  J  .......  25 

8.  Kp  as  a  Function  of  W*  .........  25 


9,10,11.  K  as  a  Function  of  L/D . 26 

.H 

12..  Kp  Design  Results.  . . .  26 

15,14,15,16.  Km  Design  Results  . . 2'f 


lv 


LIST  OF  FIGURES  (Cont'd) 

Figure  Page 

0  4 

17.  Maximum  W*T*  *  S3  a  Function  of  A*, 

AP*/T*  -  0  ............ . 27 


18.  Ejector  Performance  for  the  Primary  Nozzle 

Displaced  from  Throat  of  Secondary  Nozzle  .  .  23 

0  -4 

19.  W*T*  *  as  a  Function  jf  Primary  Nozzle 

Position  for  Constant  AP*/T*  .  28 

20.  Design  and  Operating  Characteristic  Curves.  .  3^ 

21.  Maximum  Ejector  Efficiency  as  a  Function  of 

Geometry  . . 3^ 

22.  Velocity  and  Temperature  Profiles  ......  40 

23 .  Velocity  and  Temperature  Profiles  at  a 

Low  x/D  and  w*T*0*5 . 40 

24.  Centerline  Velocity  Decay  ..........  40 

25.  Centerline  Temperature  Decay . 40 

26.  Momentum  and  Flow  Rate  Profiles . 43 


LIST  OF  TABLES 


I.  Design  Curve  Calculations  for  Sample 

Design  .  .  .  .........  33 

II.  Operating  Characteristic  Calculations  for 

Sample  Design  ........  .  .  33 


v 


REFERENCES 


1.  "Exhaust-Stack  Ejectors  for  Marine  Gas  Turbine 
Installations, "  by  A.  L.  London  and  ?.  F.  Pucci, 

TR-No,  26,  Contract  Noonr  251  Task  6,  Department  of 
Mechanical  Engineering,  Stanford  University,  1955. 

2.  "Simple  Ejector  Design  Parameters,"  by  P.  F.  Pucci, 

PhD  Thesis,  Department  of  Mecharical  Engineering, 

Stanford  University,  195^. 

5.  "Gas  Tables,"  by  J.  H.  Keenan  and  J.  Kaye,  John  Wiley 
and  Sons,  1949. 

4.  "Investigation  of  Flow  in  an  Axially  Symmetrical 
Heated  Jet  of  Air,"  by  S.  Corrisin,  NACA  ACR  3L23. 

5.  "Aerodynamics  of  Propulsion,"  by  D.  Xuchemann  and 
J.  Weber,  McGraw-Hill  Book  Co.,  1953. 

6.  "A  Note  on  the  Measurement  of  Total  Head  and  Static 
Pressure  in  a  Turbulent  Stream,"  by  S.  Goldstein,  Proc . 
Royal  Society,  Series  A,  Voi.  155,  No.  188,  pp„  570-575. 

7.  "Effect  of  Turbulence  on  Air  Flow  Measurements  Behind 
Orifice  Plates,"  by  J.  N.  Nielson,  NACA  WR  L-274. 

8.  "Aerodynamic  Measurements,"  by  R.  C.  Dean,  M.I.T.  Gas 
Turbine  Laboratory,  1953* 

9.  "The  Turbulent  Mixing  of  Co-Axial  Gas  Jets,"  by  F.  Landis 
and  A.  H.  Shapiro,  Heat  Transfer  and  Fluid  Mechanics 
Institute,  1951,  PP.  133-1^7. 

10.  "Momentum  and  Mass  Transfer  in  Co-Axial  Gas  Jets,"  by 

W.  Forstall  and  A.  H.  Shapiro,  Jour,  of  Applied  Mechanics, 
December  1950,  p.  399. 


vi 


NOMENCLATURE 


English  Letter  Symbols 

p 

A  area,  ft',  A-, ,  Ap,  A  ,  Ag,  and  Am  are  flow 

areas  defined  In  the  ejector  sketches.  Figs.  1 
and  2.  Aw  Is  the  mixing  tube  wall  area. 

C^CpiC-j  coefficients  defined  by  equation  set  [13a] 

D  diameter  of  the  mixing  tube,  ft 

f  flow  friction  factor,  dimensionless 

Ffr  wall  friction  shear  force,  # 

p 

gc  32.20  ft  lb/#  sec  ,  the  reciprocal  of  the  pro¬ 

portionality  factor  in  Newton's  Second  Law; 
force,  length,  mass,  and  time  are  the  primary 
quantities. 

h  enthalpy,  ft  #/ lb  or  Btu/lb  .  h  ,  h  ,  b„  apply 

111  O 

to  the  mixed,  primary,  and  secondary  flow  streams 
respectively. 

L  length  of  the  mixing  tube,  ft 

M  molal  mass  of  gas,  29.0  lbs/lb  mol  for  air 

p 

P  gas  pressure,  #/ft 

AP  pressure  difference,  #/ft  or  Inc! ?s  of  water 

R  the  universal  gas  constant,  15^5  ft  #/(lb  mol  R) 

T  absolute  temperature,  R 

U  local  velocity  a4:  a  point  in  the  flow  cross 

section,  ft/sec 

V  bulk  average  velocity  at  a  flow  cross  section, 

ft/sec.  Vm,  V  ,  Vs  apply  to  the  mixed,  "*•  ir.ary, 

and  secondary  flow  streams  respectively  at  the 

sections  indicated  in  Fig.  1. 

v  gas  specific  volume  (l/p),  ft'Vlb 

w  flow  rate,  lbs/sec.  wm,  wp,  ws  apply  to  the 

mixed,  primary,  and  secondary  flow  streams  re¬ 
spectively;  see  Figs.  1  3nd  2. 

x  flow  length  position  along  the  mixing  tube  mea¬ 

sured  from  the  throat  of  the  secondary  nozzle,  ft 


vii 


Greek  Letter  Symbols 


A  on  ?  denotes  a  pressure  differential 

0  denotes  a  functional  relation 

rj  energy  efficiency,  dimensionless 

p  gas  density,  lbs/ft^ 

H  gas  viscosity,  lbs/sec  ft 

Nondlmenslonal  Parameters  for  Ejector  Performance 

A*  the  area  parameter,  secondary  to  primary  flow 
area  ratio 

Km  momentum  correction  factor  relating  the  actual 

flow  stream  momentum  rate  to  that  calculated  from 
the  bulk  average  velocity,  see  Eq.  [5] 

K  static  pressure  correction  factor  relating  the 

p  actual  pressure  in  the  exit  plane  to  the  static 

pressure  at  the  ejector  wall  at  the  exit,  see  Eq.  [4] 

K*  exit  static  pressure  variation  coefficient 
P 

L/D  length  parameter  for  the  ejector 

AP*  pressure  coefficient,  a  measure  of  ejector 
performance  when  correlated  with  W* 

T*  the  temperature  parameter,  secondary  to  primary 
stream  temperature  ratio 

W*  the  flow  rate  parameter,  secondary  to  primary  stream 

flow  rate  ratio,  a  measure  of  ejector  pei  ’ormance 
when  correlated  with  AP* 

Miscellaneous 

§  denotes  pounds  force  in  distinction  to 

lbs  denoting  pounds  ma ss 
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INTRODUCTION 


The  current  status  of  ejector  design  information  repre¬ 
sents,  primarily,  the  fruits  of  a  "cut-and-try"  approach  to 
the  problem  of  ejector  development.  The  lack  of  a  systematic 
and  rational  design  method  instigated  the  present  investiga¬ 
tion,  of  which  the  preliminary  results  were  sported  in  (l) 
and  (2)*„  This  investigation  has  been  focused  on  a  specific 
application  —  the  use  of  the  exhaust  of  a  gas  turbine  as  the 
driving  medium  to  pump  a  supply  of  ventilating  air  for  the 
machinery  space.  This  application  is  characterized  by: 

a.  A  medium  subsonic  velocity,  high  temperature,  driv¬ 
ing  fluid  (exhaust  gases)  of  essentially  the  same 
molecular_welght  as  the  pumped  secondary  fluid 
(air) . 

b.  A  secondary  stream  of  ventilating  air  at  low 
temperatures  (80  to  120°P). 

c.  A  ratio  of  secondary  to  primary  flow  areas  in  the 
range  1  to  10. 

d.  A  circular  mixing  duct  of  constant  cross  sectional 
area  in  which  the  mpmentum  transfer  between  the 
primary  and  secondary  fluids  occurs. 

Tnese  are  the  characteristics  of  prime  mover  installations  in 
test  cells,  vehicles,  and  ships  incorporating  ejector-driven 
ventilating  and  cooling  systems. 

The  preliminary  results  of  the  model  tests  of  an  ejector 
with  these  characteristics  were  reported  in  (l)  and  (2).  This 
initial  study  established  the  following: 

1.  An  analytical  relation  for  the  prediction  of  ejector 
performance.  The  incoming  primary  and  secondary  velocity  and 
temperature  distributions  are  assumed  uniform  ever  the  primary 
and  secondary  areas,  respectively.  The  nonunlform  velocity 
and  temperature  profiles  at  the  exit  of  the  ejector  mixing 


* 

Numbers  in  parentheses  refer  to  the  References. 
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tube  are  accounted  for  cy  a  nondimensional  momentum  correction 
factor,  ,  modifying  the  momentum  rate  based  on  the  bulk 
average  mass  velocity. 

2.  The  relations  between  the  K  factor  and  the  design 
variables  --  area  ratio,  flow  rate  ratio,  temperature  ratio, 
and  duct  length  to  diameter  ratio.  These  relations,  established 
by  model  tests,  are  required  for  the  solution  of  the  ejector 
design  problem. 

3.  An  analysis  tentatively  indicating  that  the  number  of 

governing  nondimensional  parameters  could  be  reduced.  The 

secondary  to  primary  temperature  ratio,  T*  ,  was  combined  with 

secondary  to  primary  flow  rate  ratio,  W*  ,  in  the  form 

0.  S 

W*T*  '  to  provide  a  correlation  for  the  factor.  It  was 

also  seer  that  if  these  two  parameters  were  combined  as 
0  4 

W*T*  ’  ,  the  nondimensional  ejector  performance  parameter, 

AP*/T*  was  successfully  correlated.  In  this  way  T*  was 
eliminated  as  a  separate  parameter;  as  this  conclusion  was 
based  on  analysis,  experimental  confirmation  is  needed. 

The  present  study  Is  an  extension  of  the  foregoing  work, 

reported  in  (l)  and  (2),  with  the  following  objectives: 

0  S 

1.  Experimental  justification  for  W*T*  as  a  cor- 
; elating  parameter  for  the  momentum  correction 
factor  K_  . 

m  0  4 

2.  Experimental  justification  for  W*T*  *  as  a  cor¬ 
relating  parameter  for  ejector  performance  AP*/T*  . 

3.  A  refinement  of  the  previous  modified  one-dimensional 
analysis  with  the  introduction  of  a  static  pressure 
correction  factor  K  to  account  for  the  nonuniform 

Jr 

static  pressure  at  the  ejector  exit. 

■4.  Development  of  a  less  tedious  and  more  accurate  method 
of  determining  the  momentum  correction  factor  Km 
instead  of  the  previously  used  Pitot  tube  and  tempera- 
ture  probe  traversing  procedures  followed  by  an  in¬ 
tegration  of  the  momentum  flux  over  the  cross  section. 

5.  A  revision  and  extension  of  the  design  curves  pre¬ 
sented  in  (l)  as  a  result  of  the  investigations  under 
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3  and  4. 

Some  performance  results  for  a  combination  of  con¬ 
stant  area  and  constant  pressure  mixing. 
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SUMMARY  AND  CONCLUSIONS 


The  analysis  and  discussion  of  the  model  test  results, 
which  are  presented  In  the  later  sections  of  this  work,  lead 
to  the  following  conclusions: 

1.  It  is  empirically  demonstrated  that  the  momentum 
correction  factor  Km  (defined  by  Eq.  [5])  Is  effectively 
correlated  by  combining  the  flow  rate  and  temperature  par¬ 
ameters,  W*  and  T*  ,  in  the  form  w*T*^*^  (Figs,  6  and  7)- 

0  s 

The  correlation  of  with  W*'T*  ,  as  opposed  to  W* 

alone,  reduces  the  spread  of  the  data  by  factors  of  2  to  1 
to  10  to  1.  This  further  reinforces  the  hypothesis  regard¬ 
ing  1V*T*°  ^  as  a  correlating  factor  (l). 

0  4 

2.  It  is  empirically  demonstrated  that  W*T*  *  ef¬ 
fectively  correlates  AP*/T*  (Figs.  4  and  5);  this  verifies 
tne  conclusion  reached  analytically  that  W*  and  T*  can  be 
combined  as  W*T*0,4  for-  all  T*  >  0.3  ,  (l).  The  difference 
between  the  K  correlating  factor  (W*T*0*'1)  and  the  aP*/T* 
correlating  factor  (W*!4,  *  )  is  discussed  in  the  text. 

3.  A  static  correction  factor,  Kp  ,  which  accounts  for 
the  nonuniformity  of  static  pressure  at  the  jet  exit,  is  in¬ 
troduced  into  the  one-dimension3l  design  equation.  K  is  de- 
fined  by  Eq.  [A],  The  departure  of  Kp  from  unity  is  a  mea¬ 
sure  of  its  Influence  on  performance.  For  L/D  <  6  and 

A*  >  5  (Fig.  12),  the  effect  on  AP*  can  amount  to  20  per 
cent.  The  use  of  both  Km  and  Kp  in  the  design  equation 
allows  an  accurate  prediction  of  ejector  performance.  (In 
Figs.  4  and  5,  the  influence  Is  minor.) 

4.  A  more  accurate  method  of  determining  Km  by  ex¬ 
periment  for  low  L/D's  is  to  compute  it  from  the  performance 
results  using  the  design  equation  [13].  This  procedure  differs 
from  the  previous  one,  (2),  of  evaluating  Km  from  velocity 
and  temperature  traverses  of  the  mixing  tube  flow.  The  un¬ 
certainty  in  K  is  reduced  substantially  for  L/D's  less 

^  m 

than  4.  For  larger  L/D's  ,  the  traverse  method  yields  ac¬ 
curate  results.  The  overall  performance  method  is  used  in 
the  evaluation  of  K  for  Figs.  13,  14,  and  15. 

ITi 
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5.  On  the  basis  of  the  present  model  test  results,  the 
Km  design  curves  of  (l)  are  revised  and  extended  for  L/D's 
of  l,  2,  and  3  (Figs.  13,  14,  15).  For  L/D's  greater-  than 
4,  no  revision  was  necessary;  here  Fig.  lG  is  Fig.  21  of  (l). 
The  revision  for  L/D's  of  1,  2,  and  3  was  necessary  due  to 
a  change  in  flow  pattern  as  discussed  in  the  text.  The  K„ 
results  are  presented  as  a  function  of  A*  ,  with  W*T* as 
a  parameter.  The  range  of  A*'s  tested  extends  from  1  to  7, 
the  range  of  L/D's  1  to  8.  This  range  of  A* 1 s  and  L/D's 
will  include  most  ejector  design  applications  for  which  the 
primary  (driving)  stream  is  of  necessity  a  low  or  medium  sub¬ 
sonic  speed  flow  (e.g.  exhaust  from  a  gas  turbine,  diesel,  or 
similar  prime  mover). 

0  4 

6.  The  maximum  W*T*  *  possible  in  an  ejector  is  pre¬ 
sented  as  a  function  of  A*  ,  with  L/D  as  a  parameter  (Fig. 
17)0  This  characteristic  is  determined  from  the  design  equa¬ 
tion  [13]  with  AP*/T*  =  0  .  These  characteristics  define  the 
limits  of  ejectur  flow  rate  performance  and  aid  in  specifying 
the  conditions  for  W*  and  A?*  for  maximum  efficiency  (see 
Application  to  Design). 

7.  Some  limited  results  are  presented  for  ejector  per¬ 
formance  with  the  primary  nozzle  discharge  displaced  from  the 
throat  of  the  secondary  nozzle  (Fig,  18).  For  the  A*  tested 
(1.835),  gains  of  3$  to  10$  in  w*T*0,/;  were  realized.  The 
maximum  gains  occurred  at  the  highest  W*T*0,1+,s  .  The  prim¬ 
ary  nozzle  discharge  was  displaced  about  one  mixing  tube  di¬ 
ameter  from  the  throat  of  the  secondary  nozzle  into  the  plenum 
chamber  to  obtain  this  maximum.  The  optimum  distance  into 

the  plenum  chamber  decreased  as  w*T*0,i|  was  decreased.  The 
variation  of  w*T*^*^  with  primary  nozzle  displacement  re¬ 
vealed  that  the  ejector  performance  is  relatively  insensitive 
to  small  displacements  of  the  primary  nozzle  from  the  throat 
of  the  secondary  nozzle  (Fig.  19). 

8.  An  ejector  design  procedure  illustrating  the  use  of 
the  design  curves  is  presented  in  the  APPLICATION  TO  DESIGN 
~r  •  ■  page  29. 


5 


9.  It  is  believed  that  the  simple  ejector  problem  has 
been  satisfactorily  solved  in  its  more  essential  features, 
at  least  from  the  designers'  viewpoint.  The  evident  research 
problems  remaining  are: 

a.  Determination  of  the  optimum  primary  nozzle 
location . 

b.  Effect  of  secondary  nozzle  shape  on  ejector 
performance  together  with  primary  nozzle 
location. 

c.  Analytical  study  of  the  mixing  process,  with 

the  prediction  of  velocity,  temperature,  and 

pressure  distributions  as  a  goal  so  as  to 

supplant  the  present  empirical  procedure  of 

establishing  K  and  K 
m  p 
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A  simple  analysis  follows,  closely  paralleling  that  pre¬ 
sented  in  (i),for  the  ejector  system  described  by  Fig.  1.  The 
ejector  primary  stream  discharge  is  the  throat  location  con¬ 
figuration  of  Fig.  2a;  the  primary  nozzle  discharge  coincides 
with  the  throat  of  the  secondary  nozzle.  The  following  analy¬ 
sis  is  presented  specifically  for  this  geometry. 

The  primary  flow,  w  ,  enters  the  mixing  tube  at  section 

1  through  a  nozzle  of  area  A  with  a  velocity  V  .  This 

P  P 

primary  flow  induces  a  secondary  flow,  w„  ,  which  enters  sec- 

o 

tion  1  through  an  annulus  of  area  A  with  a  velocity  V 

s  s 

The  resulting  mixed  flow,  wm  =  w  +  wg  ,  leaves  the  mixing 
tube  at  section  2  with  mean  velocity  V  .  Tne  primary  veloc- 
ity,  V  ,  3s  greater  than  the  secondary  velocity,  V  .  The 

t'  s 

corresponding  temperatures  and  densities  are  also  shown  in 
Fig.  1. 

The  equations  for  ejector  performance- are  established  by 
combining  the  equations  of  continuity,  energy,  momentum,  ana 
state,  the  Second  Law  of  Thermodynamics,  the  defining  equations 
for  Km  and  Kp  ,  and  the  prescribed  boundary  conditions.  In 
addition  the  following  idealizations  and  definitions  are  em¬ 
ployed  : 

1.  Ste.ady-stato  flow  between  sections  0  to  1  and  1  to 

2  (Fig.  1). 

2.  Adiabatic  flow  throughout  with  one-dimensional  isen- 

troplc  flow  (reversible  adiabatic)  of  the  secondary  stream 

from  section  0  to  1,  and  irreversible  adiabatic  mixing  of  the 

primary  and  secondary  flow  streams  between  section  1  and  2. 

5.  At  section  1,  the  stream  velocities,  V  and  V„  , 

P  0 

and  temperatures,  T_  and  Tc  ,  are  uniform  across  their  re- 

p  - 

spective  streams.  However,  Yp  is  not  equal  to  Vg  and  Tp 
is  not  equal  to  T  ,  (Fig.  1 ) 

4.  Nonuniform  accelerations  and  deaccelerations  result 
in  a  nonuniform  static  pressure  distribution  across  the  stream 
at  sections  1  and  2.  This  static  pressure  distribution  is 
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accounted  for  by  a  nondlmensional  pressure  correction  factor, 

Kp  ,  relating  the  actual  static  pressure  acros."  the  mixing 
tube  to  the  wall  pressure  at  that  section.  It  was  experimen¬ 
tally  determined  that  the  wall  pressure  at  exit  is  equal  to 
atmospheric  pressure;  hence  K  F,  =  P  .  The  specification  of 
one-dimensional  lsentropic  flow  from  section  0  to  1  makes 

K  =  1  at;  section  1  since  the  static  pressure  is  uniform  in 
P 

the  plenum  chamber  (section  0). 

5.  Incomplete  mixing  of  the  primary  and  secondary 
streams  between  sections  1  and  2  and  wall  friction  resuxt 

in  a  nonuniform  velocity  distribution  which  may  be  accounted 

for  by  the  use  of  a  momentum  correction  factor,  K  , 

m 

(nondlmensional)  which  relates  the  actual  momentum  rate  at 
section  2  to  the  nominal  rate  based  on  the  bulk  average  ve¬ 
locity  and  density  at  section  2. 

6.  Both  flows  behave  as  perfect  gases  and,  moreover, 
their  molal  masses  (molecular  weights)  and  other  thermal  prop¬ 
erties  approximate  those  of  air  at  the  appropriate  temperatures. 

7.  Wall  friction  between  sections  1  and  2  is  accounted 
for  with  a  conventional  friction-factor  term. based  on  the  bulk 
flow  velocity  V  and  the  mixing  tube  wall  area  Aw  (Eq. 

16]). 

8.  Pressure  changes,  Pq  to  P^  ,  and  to  ?£  ,  are 

small  relative  to  the  static  pressure  so  that  the  gas  density 
is  essentially  dependent  on  temperature  (and  atmospheric 
pressure ) . 

9.  Flow  potential  energy  of  position  and  kinetic  energy 
changes  are  relatively  negligible  for  the  purpose  of  evaluat¬ 
ing  the  bulk  average  temperature  at  section  2. 

The  Continuity  Equation  --  The  conservation  of  mass 
principle,  applied  to  the  mixing  tube  between  sections  1  and 
2  as  a  control  volume,  yields : 


where 
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vv 


m 


=*  (p 


A  V  )  q 
m  m  m '  2 


w_ 


-  (ppVp’i 


w. 


-  (PsVs'l  ‘lal 


The  Energy  Equation  --  The  conservation  of  energy  princi¬ 
ple,  applied  to  the  nixing  tube  between  sections  1  and  2  as  a 
control  volume,  yields: 


wmhm  =  wphp 


+  wshs 


[2] 


where  the . potential  energies  of  position  and  the  kinetic  en¬ 
ergy  terms  are  neglected  (idealization  9) ,  and  the  flow  is 
adiabatic  (idealization  2). 

From  the  specification  of  perfect  gas  behavior  (idealiza¬ 
tion  6) : 

Tpl  =  0(h.n)  only  , 


from  equation  of  state  considerations.  Therefore,  with  the 
primary  and  secondary  flow  rates  and  temperatures  known,  and 
with  the  use  of  available  equation  of  state  data  (3),  Tm 
can  be  calculated. 

The  Momentum  Equation  —  The  conservation  of  momentum 
principle,  applied  to  the  mixing  tube  between  sections  1  and 
2  as  a  control  volume,  yields: 


w  V  w  V 
p  p  s 

e,  + 


P,A,  = 

11  m 


w  V 
mm 


+  (Kn  PJA, 


-4-  ? 

+  4fr 


[3] 


From  the  specified  geometry,  A^  =  A2  (Fig.  l).  Idealization 
3  allows  the  use  of  (w  V  )/g„  and  (w  V  )/g  as  the  pri- 

r.iary  and  secondary  momentum  rates  at  section  1,  respectively. 
Specifications  4  and  5  allow  the  introduction  of  Kp  and  Km 
to  express  the  pressure  force  and  momentum  rate  at  section  2. 

The  Static  Pressure  Correction  Factor  --  Kp  ,  in  Eq.  [3] > 
is  defined  by: 

where  Pwall  -  Patm  when  L/D  =  x/D  . 


h 

=  T 


wall 


/  p  <«<*-) 
Am 


m 
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The  Momentum  Correction  Factor  —  X 


defined  by: 


m 


in  Eq, 


k  - 

m 


PmV 


m  m 


/  Pu2d(|-) 

Am 


[5] 


where  Vm  is  the  bulk  average  velocity  defined  by: 


w 


m 


V  = 
m  P.JV 

1  m 


[lb] 


m  m 


The  Perfect  Gas  Equation  of  State  is  used  to  evaluate 


rr 


Pm 


(r/m)t 


[2a] 


m 


where 

[2]. 


T 


m 
KppS 


is  the  bulk  mean  temperature  evaluated  from  Eq. 

is  the  wall  pressure  measured  at  Section  2,  modi- 

to  yield  the  average  discharge  pressure. 

P 

The  wall  skin-friction  force  Ff^  is  related  to  the  bulk 
properties  using  idealization  7: 


flea  by  K 


F  =  f 
rfr  1 


.■all 


V 

m  rm 


2g, 


[6] 


A  reasonable  value  of  f  is  0.002  ,  an  average  of  a  pipe-flow 
type  friction  factor  and  a  flat-plate  type  friction  factor. 

The  value  of  f  is  not  critical;  a  variation  in  f  of  three 
to  one  does  not  affect  the  prediction  of  ejector  performance 
significantly  (l). 

The  ejector  "pumped  head"  performance  is  defined  in  terms 
of  the  plenum  chamber  pressure,  PQ  ,  and  the  discharge  pres¬ 


sure  =  ^a  *  With  the  assumption  of  isentropic  second¬ 

ary  flow  from  section  0  to  section  1  (idealization  2),  the 
energy  equation  in  differential  form  for  the  secondary  flow 
process  0  to  1  is:  2 

dh  =  "d(2g~) 

c 

The  kinetic  energy  at  section  0  (the  plenum  chamber)  is  neg¬ 
ligible.  As  the  flow  is  isentropic,  from  the  Second  Law 
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considerations : 


T  ds  =  dh  -  vdP  =  0 
and  for  small  pressure  differences 

v  4  i/p  «  i/pg 

a  constant.  Combining  these  relations  and  integrating  from 
section  0  to  1  yields: 

P  -  P.  V  2 

O  1  S  rr 


The  Solution  for  Ejector  Performance  —  The  solution 
for  ejector  performance  is  found  by  combining  the  preceding 
equations  to  yield: 

1  w2  w2  1  A  wf.  f  A 

lV  p  \  _  •  J  P  ,  s  ■> _ E  m  ir  , _ £L 

'  a~  o'  “  n-  A  I  A  n  T  A  ^  X  “  O  A  ~  A  a  t  O  S 


g  A„  A^o  A  o 
°c  m  prp  Sr 


"s  ,  £\i 

SPS  "  2  As 


wm  K  +  ££w 

Vm  L  m  2  Am 


-  V1  -  V 


where  A  and  p  apply  to  the  primary  flow  at  section  1; 

A  and  p  ,  the  secondary  flow  at  section  1;  and  K  ,  K  , 
s  s  m  p 

P.  ,  and  p  .  the  mixed  flow  at  section  2.  This  eauatlon 
m  '  r  tn  ' 

can  now  be  used  in  design  calculations  for  ejector  performance 
If  a  diffuser  having  an  efficiency  r|D  is  introduced  at 
the  end  of  the  ejector,  an  improved  plenum  chamber  vacuum  can 
be  obtained  for  a  given  secondary  flow  rate.  For  the  ejector 
with  a  diffuser  the  foregoing  analysis  may  be  extended  to 
yield  a  modification  of  Eq.  [8]: 


(Pa-PQ)  =  v  ■>" 

a  0  gcAm  |Vp 


Vs 


1 _ 

2A. 


-  TTp-  Km  +  7  J - ? 

r*i»  m  m 


+  M1  ‘  Kp) 


Diffuser 

Term 
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where  the  diffuser  efficiency  is  defined  by: 


<Pa  -  P2> 


2 'actual 


(Pa  “  ^2 ^isentropic  flow 


[10] 


For  small  pressure  differences  relative  to  the  total  pressure, 
as  considered  here,  "pressure"  and  "energy"  diffuser  efficien¬ 
cies  are  essentially  equal: 


Nondlmenslonal  Parameters  --  Several  advantages  are  real¬ 
ized  if  the  ejector  design  equation  is  presented  In  nondimen- 
sional  form.  Design  calculations  are  more  readily  accomplished 
if  nondimensional  parameters  are  used,  and  the  number  of  in¬ 
dependent  variables  Is  reduced.  Most  important  is  the  fact 
that  the  nondimensional  parameters  are  needed  to  provide  the 
criteria  of  similarity  between  model  and  prototype  in  order 

to  establish  K  and  K  from  the  model  tests, 
p  m 

The  nondimensional  parameters  useful  for  this  analysis 


which  compares  the  "pumped  head", 

(Pa  -  p0)/ps  >  for  the  secondary  flow  to 

the  "driving  head",  V?/2g  ,  of  the 

P  c 

primary  flow. 

W*  =  ws/wp  a  flow  rate  ratio,  secondary-to¬ 
primary  mass  flow  rate. 

T*  =  T  /T  an  absolute  temperature  ratio, 
s  p 

T*  =  Pp/ps  >  a  flow  density  ratio,  since 
the  exhaust  gas  is  presumed  to  have  close 
to  the  same  molal  mass  as  the  ventilating 
air,  and  pressure  differences  are'  small 

relative  to  (?t"»osr'her1c  pressure.. 
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A  /A  ,  A  /An,  A*  £  A  /A  ;  all  flow  area  ratios  defined  by 
p  m  m  s  p 

outer  duct  inside  dimensions,  primary 
tube  ID,  primary  tube  wall  thickness, 
and  diffuser  discharge  area. 

Aw/Am  =  wall  friction  area  to  mixing  tube  flow 
cross  section  area  ratio. 

K  ,  K  ,  f,  ti_  =  the  pressure  correction  factor,  the  mo- 
p  m  'D 

mentum  correction  factor,  the  wall  fric¬ 
tion  factor,  and  the  diffuser  efficiency. 
These  design  parameters  have  been  pre¬ 
viously  considered. 

a  convenient 

nondimensional  grouping  for  the  effect  of 
nonuniform  static  pressure  across  the 
mixing  tube  exit  on  ejector  performance. 

Introducing  these  nondimensional  parameters  into  Eq. 

[93* 


For  a  fixed  ejector  geometry,  Eq.  [11]  reduces  to  the  form: 


AP*  K* 

-  =  C,  -  C0W*2T*  -  C,W*  (T*  +  1)  +  -2- 

T*  1  d  '  T* 


[12] 


where  the  C  coefficients  involve  only  geometry,  K  m  ,  f  , 
and  .  These  are  defined  in  the  equation  set  [13a].  From 
Eq.  [11]: 

=  0[W*,  T»,  A*, A  /Am,  Vv  W  V  V  ’’d1 


13 


With  the  geometry  established,  aP*/T*  becomes  a  function  of 
only  five  parameters  (W*,  T*,  Km,  f,  K^).  However,  calcula¬ 
tions  have  shown  that  a  fixed  v.lue  of  f  of  0.003  is  ade¬ 
quate  for  design  purposes.  Also,  experiments  indicate  that 
Kp  ,  for  a  given  geometry  (A*,  L/D),  can  be  treated  as  a  con¬ 
stant  (Fig.  8).  This  reduces  AP*/T*  to  a  function  of  only 
three  parameters,  W*,  T*,  and  Km  .  Moreover,  calculations 

performed  with  Eq.  [12]  show  the  W*  and  T*  can  be  combined 

0 

as  the. single  term  W*T*  *  '  for  all  T*  >  0.3  with  no  sig¬ 
nificant  error  [Fig.  5  and  Table  2,  (l)].  This  combination 
has  also  been  verified  by  experiment  (Figs.  4  and  5).  These 
considerations  reduce  the  description  of  the  performance  of 
a  particular  ejector  to  a  relation  between  three  variables. 


^  =  0[W*T*°-\  KJ 


with  a  consequent  simplification  of  the  ejector  model  test 
program  and  also  the  design  procedure. 

The  ejector  performance  Eq.  [12],  using  the  w*T*0-/| 
combination,  becomes  in  approximate  form: 


AP* 


-.c1  -  c2(w*t*0,4)2 


(W*T«°*‘4 


[13] 


and  this  equation  is  valid  for  T*  >  0.3  .  Note  that  for. 
T*  =  1  ,  Eqs.  [i2]  and  [13j  become  identical.  The  coeffi¬ 
cients  ,  and  C.^  are  defined  for  both  Eqs.  [12] 

and  [13]  as  follows: 

1  -  ^ 

m 


with 


d1  (B) 

L"lil 

2 

O) 


Km  + 
m 


2  A 


m 


2 


1  /  1  \ 
l 1  "  / 


y  [15a] 


14 


r:tal  results: 


From  the  exper1- 

Km  =  0(W *,  T*,  L/Dj  A*) 

As  it  was  expected  that  the  ratio  of  momentum  fluxes, 

KVi  _  K\2  Pg  .  ^  _  <W*T,>0~5)2 
("pVl  \Kp  j  Ps  As  A* 

was  the  controlling  mixing  parameter,  it  was  hypothesized 
that  W*  and  T*  could  he  combined  as  v.’*T*0*^  to  corre- 
late  Km  .  The  work  herein  presents  experimental  results 
that  substantiate  this  hypothesis(Flgs.  6  and  7).  Hence. 

Km  =  0(W*T*0*5,  L/D,  A*) 

and  with  the  geometry  established: 

Km  =  0(W*T*0*5) 

Experimentally,  it  was  found  that  K  varies  significantly 

A  C  .  *  *  ^ 

with  tf*T*  0  only  ior-  L/D  <•  4  .  For  L/D  >  4  ,  de¬ 
pends  essentially  on  geometry  only  (Figs.  1?,  14,  i,;.  and 
16). 

It  is  coincidental  that  the  flow  rate-temperature  cor- 

0  *5 

relation  parameter  for  ,  W*T*  ,  is  nearly  the  same  a s 

the  AP*/T*  correlation  parameter,  as  these  two 

parameters  originated  from  two  unrelated  concepts.  The  Km 

correlation  factor  was  hypothesized  on  the  consideration  that 

the  ratio  of  the  secondary  to  primary  momentum  flux  was  the 

fluid  parameter'  characterizing  the  momentum  transport  process 

0  h 

in  the  mixing  tube.  W*T*  * ^  was  accepted  as  the  correlation 
parameter  on  the  basis  of  the  model  test  results.  The  choice 
of  w*T*0,'!+  as  the  £P*/T*  correlation  parameter  resulted 
from  calculations  performed  with  the  design  equation  [12]. 

It  was  seen  that  for  T*  >  0.3  ,  the  flow  rate  and  tempera- 
ture  ratios  'could  be  combined  as  w*T*u,H  with  no  significant 
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error  in  AP*/T*  [see  Fig.  5  and  Table  2,  ( 1 ) ] .  This  al¬ 
lowed  the  replacement  of  W*^T*  by  (W*T*^ and 
V*(l  +  T*)  by  2W*T*^‘^  as  in  Eq.  [13].  Experiments  veri¬ 
fied  this  combination  as  a  correlation  parameter.  Figs.  4  and 
3.  These  combinations  of  flow  rate  and  temperature  ratios 
serve  the  important  functions  of  simplifying  the  model  test 
program,  the  presentation  of  the  .results,  and  the  design 
procedure . 


Ejector  Energy  Efficiency  --  A  thermodynamically  mean¬ 
ingful  definition  of  ejector  efficiency  is 

•  [ws(Pa  '  P0)/PSJ 


'iE 


t»p  V26j 


[«] 


‘This  is  an  expression  for  the  efficacy  of  conversion  of  the 

primary  jet  kinetic  energy  tc  secondary  stream  flow  work. 

• 

This  definition  is  analogous  to  the  usual  pump  efficiency 
comparing  the  fluid  flow  work  increase  to  the  power  input. 

In  terms  of  the  nondimensional  parameters,  [14]  becomes 

tje  =  (w*T*0*4)  (AP*/T*)T*0,6  [15] 


The  maximum  value  of  ejector  efficiency  is  useful  for  the  de¬ 
sign  problem  (see  Application  to  Design).  The  exact  expres¬ 
sion  for  maximum  ejector  efficiency  can  be  found  by  substi¬ 
tuting  the  expression  for  (aP*/T*)  in  terms *of  W*T*0’4, 

Eq.  [13],  into  [15],  differentiating  t|„  with  respect  to 

oh  o  4  ^ 

W*T*  *  ,  evaluating  the  W*T*  necessary  for  maximum  ef¬ 

ficiency,  and  inserting  this  expression  back  into  [15].  How¬ 
ever,  this  is  a  cumbersome  process,  and  unnecessary  as  a  good 
approximation  for  maximum  efficiency  can  be  readily  estab¬ 
lished  as  follows.  From  Figs.  3,  4,  and  5,  and  Fig.  16  of 

(l),  w*T*0,4  is  seen  to  be  essentially  a  linear  function 

0  4 

of  AP*/T*  .  Hence,  tje  will  be  a  maximum  when  W*T* 
equals  (W*T*0*4)  Q  /2  and  AP*/T*  equals  (AP*/T*)ma  /2  . 

ukaX  liid  A 

Also,  from  [13], 


(AP*/T*)[R0X  ■=  0,  +  KJ/T* 


£  16  j 


Therefore 


y 


(O  =  i  (W*T*°‘4)  (C,  +  k*/T*)T*0*6 

v  'E'max  4  v  'max  1  P 

This  is  the  approximate  expression  for  maximum  ejector 
efficiency. 


[IT] 


17 


FIGURE  1 

SIMPLE  EJECTOR  WITHOUT  A  DIFFUSER 

Illustrates  idealizations  regarding  the  entering 
velocity  and  temperature  profiles. 

FIGURE  2 

ALTERNATIVE  LOCATIONS  FOR  THE  PRIMARY  JET  DISCHARGE 

Location  (c)  or  (d)  gives  the  best  performance  and 
location  (b)  gives  poorer  performance.  The  analysis 
presented  in  this  work  is  for  throat  location  (a). 
Limited  test  results  are  reported  for  locations  (b) 
and  (c ) . 


FIGURE  3 

EJECTOR  PERFORMANCE  -  INFLUENCE  OF  DEPARTURES  FROM 
SIMPLE  ONE-DIMENSIONAL  THEORY 

1)  Km  =  1  ,  K  =  1  ;  simple  one-dimensional  theory. 

2)  K  >  1  ,  K  =  1  ;  one-dimensional  theory  modified 
for  nonuniform  momentum  rate  profile  leaving  the 
ejector . 

3)  K  >  1  ,  K  <  1  :  o^e-dimensional  theory  modified 
'  m  ’  P 

for  both  nonuniform  momentum  rate  and  static 
pressure  profiles  leaving  the  ejector. 


FIGURES  4  and  3 

EJECTOR  PERFORMANCE  -  COMPARISON  OF  PREDICTED  AND 

MODEL  TEST  RESULTS 

Shaded  area  indicates  uncertainty  in  predicted  results 
due  only  to  uncertainty  in  evaluation.  Note  also 

that  additional  uncertainties  exist  in  the  measurements 
of  AP*  ,  T*  ,  and  W*  .  See  Appendix  B  for  a  summary 
of  the  uncertainties. 


IS 
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MODEL  TEST  RESULTS 


K  Correlation  --  One  of  the  primary  objectives  of  this 
Investigation  was  to  determine  if  w*T*J*-:>  is  an  adequate 
correlating  factor  for  the  momentum  correction  factor  Km  . 
This  choice  of  W*T*^ ' 3  resulted  from  the  intuitive  conclu¬ 
sion  that  the  ratio  or  secondary  to  primary  momentum  flux 
characterized  the  momentum  transport  process: 


KVS>1  _  4  P£  ^  =  <W*T* 

(WPVP>1  Wp  AS 


0.5' 


A* 


To  test  this  choice,  the  expression  for  Km  at  x/D  =  0 
and  L/D  >  0  ,  designated  Km  ,  was  formulated  theoreti¬ 
cally,  (l),  with  the  idealizations  that  V,  ,  ,  T„  ,  and 

T  are  constant  across  the  secondary  and  primary  streams, 
respectively  (idealization  5,  page  7  ): 


K 


m,  o 


_ ]____  _ 1 _ 

(ApAm)  ( W*T*  +  1)(W*  +  1) 


1  + 


W*  ^  T* 

A* 


[18] 


The  specification  that  L/D  is  greater  than  zero  in  evaluat¬ 
ing  K  is -necessary  to  Inst  re  that  there  will  be  an  in- 
m,o 

duced  flow.  Here,  if  (W*T*  +  l)(W*  +  l)  is  replaced  by  the 

approximation  (l  +  W*T*  °)  ,  there  is  no  significant  error 

introduced  in  Km  q  providing  T*  >  0.5  ,  Fig.  25,  p.  5^> 

(l).  Thus  with  the  ejector  geometry  specified,  K  becomes, 

n  c; 

to  a  good  approximation,  a  function  of  W*T*  •  ^  only,  in¬ 
stead  of  W*  and  T*  separately: 


K. 


(VAm) 


(1  +  W*T*0,5)2 


!  +  (W»T»0,5)2 
A* 


[19] 


0  5 

Another  factor  influencing  the  choice  of  W*T*  is 
that  the  dependency  of  Km  on  W*  and  T*  decreases  as  the 
length  (L/D)  of  the  ejector  is  increased.  The  maximum  depend¬ 
ency  of  K  on  W*  ana  T*  occurs  at  an  x/D  of  zero,  and 

in 

for  x/D  >  4  ,  K  is  essentially  a  function  of  geometry 

m 
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0  5 

alone.  This  fact  makes  the  acceptance  of  W*T*  *  as  the 
combined  flow-temperature  parameter  an  even  better  approxi- 
mati  r.  than  suggested  above. 

To  test  these  considerations,  K  was  evaluated  at 

'  m 

x/D's  of  0  and  2  with  A**  a  of  1.835  and  6.984  (L/D  =  8) 

by  integration  of  the  velocity  and  temperature  profiles  at 

these  stations  in  the  mixing  tube.  The  resulting  K  '  s  were 

0  c.  11  ^ 

plotted  as  functions  of  both  W*  and  W*T*  *-  (Figs.  6  and 
7).  The  tests  were  carried  out  for  temperature  ratios  of 
0.74  and  1.00  .  For  the  x/D  -  0  results,  the  relation  be¬ 
tween  Km  and  W*T*^*-’  is  Eq.  [19]  for  T*  >  0 . 3  .  The 

standard  deviations  of  the  data  about  K  for  K  as  a 

m,  o  m 

function  of  V*T*  and  for  K_.  as  a  function  of  W*  were 

m 

determined.  As  no  analytical  prediction  exists  for  Km  at 

x/D  =  2  ,  least-souare  straight  lines  were  drawn  through  the 

OS 

data  for  K  as  a  function  of  W*T*  and  for  K  as  a 
m  m 

function  of  W*  .  (The  experimental  results  indicate  that 
Km  is  nearly  linear  with  W*T*0*^  for  x/D  >  2  .)  The 
standard  deviations  of  the  data  about  these  lines  were  also 
determined.  The  resulting  values  of  the  relative  standard 
deviations  (aw*/aw*T*0. 5]  are: 

A*  =  1.835  A*  =  6.984 


x/D  =0  2.1 

x/D  =  2  10.7 


5,6 


5.2 


As  the  correlation  of  K  with  w**p*°*5  reduces  the 
spread  of  the  data  by  factors  of  2:1  to  10:1,  w*T*J,!3  is 

definitely  acceptable.  However,  additional  confirmation  at 
test  magnitudes  down  to  T*  —  0.3  would  be  desirable. 


AP*/T*  Correlation  --  The  flow  rate -temperature  corre¬ 
lating  parameter  for  AP*/T*  is  W*T*0’^  .  This  choice  was 
based  on  calculations  with  the  ejector  design  equation  [12] 
which  demonstrated  that  W*  and  T*  could  be  combined  as 
for  all  T*  >  0.3  with  no  significant  change  in 
the  value  of  aP*/T*  [(l).  Fig.  5  and  Table  4],  The  results 
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no  del  tests 


of  the  model  tests  substantiate  this  conclusion  (Figs,  4  and 

5).  The  values  of  AP*/T*  for  T*  =  0.74  are  consistently 

0  4 

in  agreement  with  those  for  T*  =  1.00  :  thus  W*T*  “  is 
verified  as  the  AP*/T*  correlation  parameter. 


Design  Results  --  In  (l)  and  (2)  and  in  the  early 

m 

stages  of  this  work,  EKe  momentum  correction  factor,  K  ,  was 
evaluated  by  integration  of  the  velocity  and  temperature  pro¬ 
files  across  the  mixing  tube.  The  profiles  were  determined  by 
Pitot  tube  velocity  and  thermocouple  temperature  traverses. 

As  a  check  on  the  accuracy  of  the  traverse  results,  the  metered 
incoming  primary  and  secondary  flow  rates  were  compared  to 
these  obtained  by  the  integration  of  the  profiles.  The  results 
obtained  bj  integration  were  consistently  higher  than  the  me¬ 
tered  results.  This  error  may  well  be  attributed  to  the  high 
turbulent  velocity  fluctuations  in  the  mixing  regions;  Corrisin 
(4)  found  velocity  fluctuation  values  of  u'/Ulor>a-,  to 

0.5  in  the  mixing  region  of  a  free  jet.  The  velocity  fluctua¬ 
tions  simulate  a  more  even  distribution  of  momentum  flux  than 
actually  exists  in  the  mixing  region  (see  Appendix  B  for  an 
analysis  of  the  effect  of  the  turbulent  velocity  fluctuations 
on  K  ).  As  a  result,  the  K  values  obtained  from  the  ve¬ 
locity  and  temperature  profiles  are  lower  than  the  actual  Km's. 

The  K  uncertainty  dus  to  this  effect  and  instrument 
m 

uncertainty  is  estimated  to  be  only  it  2 %  for  x/D  >  4  ,  but  for 
x/D  <  4  ,  K  is  estimated  to  be  Ito  5#  too  low.  As  a  con- 
sequence  the  resulting  uncertainties  in  the  predicted  values 
of  AP*/T*  are  2  to  15  per  cent  (see  Figs.  4  and  5). 

The  uncertainty  in  AP*/T*  is  ±  yf>  when  AP*/T*  is  mea¬ 
sured  directly.  As  this  is  generally  iess  than  the  uncertainty 
in  AP*/T*  resulting  from  the  Km  traverse  results,  it  was 

concluded  that  K  calculated  directly  from  the  measured 

m 

AP*./T*  using  Eq.  [11]  would  provide  the  best  Km  data  for 
design  purposes.  The  uncertainty  of  Km  by  this  procedure 
is  estimated  to  be  about  ±  1$. 

Another  reason  for  calculating  Km  from  Eq.  [11] 
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is  that  the  length  of  the  mixing  tube  influences  the  veloc¬ 
ity  and  temperature  profiles  at  positions  along  the  mixing 
tube.  In  (l),  all  Km  values  we^e  computed  by  integration 
of  the  velocity  and  temperature  profiles  at  stations  along  a 

mixing  tube  of  length  L/D  -  5  .  The  X  's  obtained  at  po¬ 
rn 

sitions  along  a  long  mixing  tube  are  higher  than  the  Km’s 

determined  at  the  discharge  of  a  short  mixing  tube  (see  Fig. 

7  for  a  comparison  of  K  at  x/D  -  2  evaluated  from  trav- 

m 

erses  for  L/D's  of  8  and  2  respectively).  This  difference 
is  attributed  to  a  difference  in  flow  patterns  --  separation 
is  present  at  x/D  =2  in  the  long  mixing  tube,  but  is  not 
when  the  mixing  tube  length  is  L/D  =  2  .  These  difficulties 
are  circumvented  when  Is  evaluated  from  the  design  equa¬ 

tion  and  the  experi mentally  determined  AP*/‘f*  , 

The  K  results  presented  in  Figs.  13,  14,  and  15  have 
been  derived  directly  from  AP*/T*  measurements  as  described 
above.  The  L/D  >  4  results  are  unaffected  by  the  afore¬ 
mentioned  changes;  consequently  Fig.  16  is  taken  directly  from 
(1). 

These  Km  results  are  improved  over  those  presented  pre¬ 
viously  (1)  in  that:  a)  the  Km  values  have  been  improved  by 
use  of  the  new  computing  technique,  and  b)  the  results  for 
L./D  =  1,  2,  an.  p  have  been  extended  from  an  A*  range  of 
1  to  4  to  a  range  of  1  to  7. 

Kp  Design  Results  --  The  factor  is  needed  to  account 

for  the  conversion  of  kinetic  energy  to  flow  work  that  occurs 
after  the  ejector  has  exhausted  to  the  atmosphere.  Even  though 
this  process  occurs  outside  the  ejector,  it  affects  the  per¬ 
formance  in  that  the  mean  exit  static  pressure  of  the  ejector 
is  slightly  less  than  atmospheric.  This  diffusing  action  is 
generally  small  but  produces  up  to  a  20  per  cent  increase  in 
performance  in  AP*/T*  at  high  mignltudes  of  and 

also  in  W*T*^’Z|  at  high  magnitudes  of  AP*/T*  „  These  ef¬ 
fects  are  evident  in  Fig.  3* 

The  variation  of  X  with  W*  is  plotted  in  Fig.  8  for 

P 
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A*  =  1.835  •  L/D  was  taken  as  the  dependent  variable  and 

W*  was  omitted  as  the  variation  in  Kp  is  slight,  and  as 

AP*/T*  does  not  depend  markedly  on  Kp  .  It  is  seen  from 

Eq.  [Ip]  that  a  change  of  20  per  cent  in  (l  -  K  )  produces 

P 

only  a  2  per  cent  change  In  •  The  variation  or  Kp 

with  L/D  for  A* 1 s  of  1.833,  3.741,  and  6.984  is  presented 
in  Figs.  9,  10,  and  11.  Figure  12  is  the  cross-plotted  re¬ 
lation  between  K  and  A*  providing  the  design  curve  to 

]? 

accompany  Figs.  13  to  16  for  Km  . 


Nozzle  Location  Results  --  Some  limited  results  are  pre¬ 
sented  for  the  ejector  configuration  with  the  primary  nozzle 
discharge  aisplaced  from  the  throat  of  the  secondary  nozzle 
(see  Fig.  2c).  The  one  geometry  tested  was  A*  =  i.835  and 
L/D  =  8  .  Figure  18  presents  the  variation  of  with 

AP*/T*  .  At  a  given  AP*/T*  the  flow  rate  is  increased  about 
10  per  cent  at  the  low  values  of  aP*/T*  and  2  per  cent  at 
the  high  values.  The  maximum  increases  occur  when  the  nozzle 
is  one  mixing  tube  diameter  into  the  plenum  chamber 
[ (x/D)n  -  -1]  at  low  AP*/T*  and  at  (x/D)n  =  -0.25  at  high 
AP*/T*'s  .  The  variation  of  w*T*°‘^  with  (x/D)n  for  con¬ 
stant  AP*/T*  Is  presented  in  Fig.  19*  The  curves  are  quite 
flat,  indicating  that  the  position  of  the  primary  nozzle  is 

not  too  critical.  The  upstream  position  for  optimum  perform- 

0  4 

ance  is  seen  to  decrease  as  W*T*  ’  decreases. 

This  behavior  may  be  explained  by  the  fact  that  the  angle 

at  which  the  jet  spreads  is  proportional  to  (Vp  -  Vs)/(Vp  +  Vc 

(5).  Decreasing  W*  increases  (V  -  V  )/(V  +  V  );  the  jet 

p  s  p  s 

spreads  faster  at  low  W*  .  If  the  primary  nozzle  is  too  far 
into  the  plenum  chamber,  the  primary  jet  will  be  diverted  in¬ 
to  the  plenum  chamber  by  the  secondary  nozzle  walls. 

0  4 

The  fact  that  an  improvement  in  W*T*  *  occurs  by  dis¬ 
placing  the  primary  nozzle  discharge  into  the  plenum  chamber 
can  be  rationalized  with  the  help  of  the  momentum  relation. 

With  the  primary  nozzle  in  the  plenum  chamber,  the  incoming 

momentum  rate  is  K  .  [ (W  V  )/g„ 1  instead  of 

m,  1  mm  “c 
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[(wpvp)/gc]  +  [(wsVs)/gcJ.  Equation  [5]  becomes 


K 
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Now,  as  some  mixing  has  taken  place  between  the  nozzle  dis¬ 
charge  and  section  1,  Km  ^  is  less  than  Km  ;  Eq.  [5] 

could  oe  written  with  K  [ (w  V  )/g  ]  instead  of 

m,  0  mm  c 

[ (w  V  )/g  ]  +  [(w  V  )/g  ]  .  Hence,  for  the  same  AP*/T*  , 
ppc  ssc 

which  is  proportional  to  (K^P^  “  Pi  )  1  W*T*  *  must  increase. 
While  this  effect  should  be  investigated  further,  a  tentative 
conclusion  is  that  a  performance  gain  of  about  10  per  cent,  may 
be  obtained  by  moving  the  primary  nozzle  about  one  mixing  tube 
diameter  into  the  plenum  chamber. 

An  investigation  of  the  influence  of  the  shape  of  the 
secondary  flow  nczzle  together  with  primary  nozzle  location 
may  demonstrate  that  even  greater  improvements  in  performance 
are  possible,  (see  Fig.  2d). 
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FIGURE 


r 


K  AS  A  FUNCTION  OF  W*  AND  W^T*0*5 
m 

Demonstrates  correlation  provided  by  combining  W* 

and  T*  in  the  form  K  Is  derived 

m .  o 

theoretically  from  the  idealizations  of  uniform  flow 
velocity  and  temperature  of  both  the  primary  and 
secondary  streams  at  Section  1. 


FIGURE  7 

K  AS  A  FUNCTION  OF  W*T*°“5 
m 

Demonstrates  correlation  provided  by  combining  W* 
and  T*  in  the  form  W*T*0,'i.  Also  illustrates  the 
influence  of  overall  length  of  the  mixing  tube  at 
x/D  ==  2  ,  where  L/D  =  2  is  compared  to  L/D  =  8  „ 
The  x/D  =  0  line  is  the  analytical  derived  relation 
equation  [19].  The  x/D  a  2  line  is  based  on  ex¬ 
periments  with  T*  =  1  . 


FIGURE  8 

Kp  ASA  FUNCTION  OF  W* 

Illustrates  the  influence  of  L/D  .  The  results 
are  obtained  from  integration  of  the  static  pressure 
profiles . 
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FIGURES  9,  10,  and  11 


K  AS  A  FUNCTION  OF  L/D 

r-' 

See  Page  23  for  the  effect  of  using  the  average 
value  of  in  the  design  equation  [13], 


FIGURE  12 


Kp  DESIGN  RESULTS 


Curves  are  cross-plotted  from  Figs.  9,  10,  and  11 


L/D 


FIGURES  13,  14 >  13,  and  16 

K  '  DESIGN  RESULTS 
m 

Figures  13,  1-4,  and  15  are  to  be  used  in  place  of 
Figs.  18,  19,  and  20  of  TR-No.  26  (l).  Figure  16 
is  Fig.  21  of  TR-No.  26. 


FIGURE  17 

MAXIMUM  W*T*0*^  AS -A  FUNCTION  OF  A*  ,  AP*/T*  =  0 

The  values  are  computed  from  the  design  equation 
[13]  and  from  Figs.  13,  14,  15  and  16. 
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*1  *)ls 


FIGURE  18 


EJECTOR  PERFORMANCE  FOR  THE  PRIMARY  NOZZLE  DISPIACED 
FROM  THE  THROAT  OF  THE  SECONDARY  NOZZLE 


FIGURE  19 

T*0,4  AS  A  FUNCTION  OF  PRIMARY  NOZZLE  POSITION,  (x/D)n 

FOR  CONSTANT  AP*/T* 


APPLICATION  TO  DESIGN 


The  information  presented  in  this  work  enables  the  de¬ 
signer  to  define  the  geometry  of  a  simple  jet  ejector  as¬ 
suring  delivery  of  a  specified  secondary  flew  under  specified 
operating  conditions.  The  design  is  formulated  as  follows: 

1.  The  design  point  specif ica cions  of  the  ejector  are 
first  listed.  These  include  the  primary  flow  temperature, 
mass  flow  rate,  flow  area,  density,  and  exit  pressure.  These 
quantities  will  be  fixed  by  the  source  of  the  primary  stream, 
probably  the  exhaust  of  a  prime  mover.  Additionally,  the  de¬ 
sired  secor.dary  air  flow  mass  rate  and  temperature  will  have 
to  be  established  from  the  cooling  or  ventilating  require¬ 
ments.  An  estimate  of  the  pressure  rise  for  the  ejector  is 
also  required;  this  estimate  is  determined  by  the  ejector 
exit  pressure  relative  to  the  minimum  plenum  chamber  pressure. 

2.  From  the  values  of  the  flow  rare  and  temperature 

specif ications,  w  ,  w  ,  T  ,  and  T  ,  and  the  estimated 
pressure  rise,  listed  under  step  1,  the  nondimensional  per¬ 
formance  parameters,  (W*T*®‘^),  „  and  (AP*/T*)  ,  at 

a.p.  a.p. 

the  design  point  are  computed. 

5,  The  value  of  (W*T*0'"*)  is  now  computed  by  doub- 

r\  2i  ITi3X 

ling  (W*T*  ),  (see  Analysis).  Figure  17  is  entered 

Cl  *  p  # 

with  this  value,  and  the  compatible  geometries  (L/D  and  A*) 
are  listed.  The  values  of  (A?*/T*)  v  corresponding  to 
these  geometries  are  evaluated  using  Eq.  [16];  is  de¬ 

fined  by  tne  equation  set  [15a},  the  values  of  Km  are  taken 
from  the  appropriate  Figs.  15,  14,  15  or  16,  and  the  values 

of  K  from  Fig.  12.  The  values  of  (AP*/T*)  and  the 
P  0  4,  mdX 

value  of  (W*T*  *  ;  are  entered  on  a  design  plot  of 

0  4  max  ✓ 

W*T*  "  as  a  function  of  AP*/T*  ,  as  in  Fig.  20.  The  de¬ 
sign  characteristics  for  the  various  geometries  are  then  ap- 

0  4 

proximate!  as  linear  relationships  between  W*T*  '  and 
AP*/T*  .  The  design  point  determined  in  step  2  is  entered 
on  the  design  plot.  The  geometry  of  the  ejector  'that  will 
yield  the  desired  performance  is  defined  by  the  design 
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characteristic  that  passes  through  the  design  point.  (This 
procedure  will  yield  a  design  with  a  slightly  higher  value 
of  AP*/T*  than  the  design  value  due  to  the  actual  slight 
convex  up  curvature  of  the  operating  characteristic.) 

4.  With  the  geometry  determined,  a  more  accurate  oper¬ 
ating  characteristic  is  established  with  the  aid  of  Eqs.  [13] 
and  [13a],  Fig.  12,  and  the  pertinent  Figs.  13,  14,  15,  or 
16.  The  operating  point,  determined  by  the  intersection  of 
the  line' of  constant  (W*T*^*^)^  and  the  operating  char¬ 
acteristic,  is  noted  on  the  design  plot,  and  the  correspond¬ 
ing  values  of  (AP*/T*)  and  are  computed. 

O  .  p  •  Hj 

In  order  to  illustrate  the  above  design-  method,  a  sample 
design  will  be  carried  out  for  the  gas  turbine  installation 
described  in  Table  1,  (l). 

1.  Primary  and  secondary  flow  specifications: 

Primary  flow  rate  w  7. 00  lb/sec 

v 

Primary  flow  temperature  Tp  =  1660  R 

Exhaust  back  pressure  P  =14.7  psia  =  1  atm 

P  0 

Primary  flow  area  Ap  =  1.07  ft 

Primary  flow  density  PD  =  0.0239  lb/ft'' 

Primary  flew  velocity  Vp  =  274  ft/sec 

Secondary  flow  rate  w  =  3*50  Ib/sec 

o 

Secondary  flow  temperature  T  =  390  R 

S  p 

Secondary  flow  density  ps  =  O.O672  lb/ft 

Desired  minimum  ejector  0 

pressure  rise  AP  =  l'^O  =  3.00  #/f t^ 


Design  point  values ; 


0.330 


30 


3-  The  geometries  that  will  deliver  the  required 
(Vs*T*°‘  ')  -  2(W*T*0,tt)  =  C.660  ,  and  the  values  of 

IlUi  X  d  «  jj  • 

(AP*/T»)  associated  with  these  geometries,  are  presented 
In  Table  I.  The  design  curves  corresponding  to  these  geometries 
are  plotted  on  Fig.  20.  The  iesign  curve  that  falls  closest 
to  the  design  point  is  curve  b);  the  geometrical  constants 
are  L/D  =  3  and  A*  =  2.31  .  The  operating  characteristics 
can  now  be  determined  more  accurately  using  Ea.[13];as  K  is  a 
function  of  W*T*U°  Fig.l  is  also  used. 

4.  Tne  data  needed  to  establish  the  operating  charac¬ 
teristic  is  presented  in  Table  II The  resulting  ejector  op¬ 
erating  characteristic  is  shown  on  the  design  plot.  It  is 
seen  that  the  operating  point  occurs  at  a  slightly  higher 
AP*/T*  than  the  design  point  value.  The  operating  point 
values  are: 

(W*T*0*4)OjP>  =  0.330 
( AP*/T* )  =  0.197 

O  o  P  • 

The  ejector  operating  efficiency,  computed  from  Eq.  [15],  is: 

Hg  =  3.5  per  cent 

The  above  design  procedure  defines  the  ejector  of  minimum 
length  (L/D  =  3  ,  A*  =  2.31)  that  will  provide  the  desired 
performance.  If  a  longer  ejector  Is  allowable,  and  a  conser¬ 
vative  design  is  desired,  a  design  curve  such  as  c )  of  Fig. 

20  may  be  chosen  (L/D  =  4  ,  A*  =  1.88).  The  operating  condi- 

0  4 

tion  will  now  deliver  che  same  W*T*  ’  at  a  higher  value  of 

AP*/T* ,  alternatively,  a  higher  W*T*  at  the  design  point 

value  of  A?*/T*  may  be  obtained. 

As  previously  mentioned,  the  concept  of  ejector  efficiency 

is  useful  and  thermodynamically  meaningful  In  that  It  provides 

a  basi~  of  comparison  between  the  ejector  and  the  conventional 

fluid  pump.  The  ejector  definitely  suffers  eff iciency-wise 

/  0  6 

in  comoarison,  as  can  be  seen  in  Fig.  21  where  /T*  * 

’  °  ‘E  max 

as  a  function  of  A*  for  various  L/D's  is  presented.  The 
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maximum  efficiency  realizable  for  the  ejector  characteristics 
presented  in  this  work  is  about  12  per  cent.  While  the  ef¬ 
ficiency  is  low,  it  is  acceptable  if  the  ejector  exploits  en¬ 
ergy  that  would  otn^rwise  be  wasted.  The  ejector  eliminates 
the  initial  cost,  weight,  and  power  requirements  of  the  blow¬ 
ers  that  would  be  used.  Also,  an  increase  in  plant  thermal 
efficiency  is  realized  by  the  reduction  in  exhaust  back  pres¬ 
sure.  For  a  gas  turbine  plant,  these  considerations  can  re¬ 
sult  in  a  net  saving  of  one  to  one  and  one-half  per  cent  on 
fuel  consumption,  (l). 
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TABLE  I 


Design  curve  calculations  for  sample  design 


Curve 

L/D 

A  * 

K 

m. 

C1 

Ko 

K*/T* 

P 

(AP*/T*) 

'  '  'max 

a) 

2 

2.97 

2.08* 

0.238 

0.99945 

0.041 

0.279 

b) 

- 2 i 

y 

2.31 

1.37* 

0.350 

0.99970 

0.022 

0.372 

c) 

4 

1.88 

1.15 

0.412 

0.99975 

0.019 

0.431 

d) 

6 

1.69 

1.03 

0.447 

0.99990 

0.007 

0.454 

e) 

8 

1.53 

1.015 

0.459 

0.99995 

0.004 

0.463 

*As  Km  for  W*T*0,ij  =  0  is  not 
evaluated  for  these  geometries,  a 
value  for  Km  slightly  higher  than 
that  for  =  0.2  is  chosen 

from  Figs.  14  and  15. 


TABLE  II 


Operating  characteristic  calculations  for  sample  design 


W*T*^ 0  ^ 

. 5 

K 

m 

6 

Ci 

C2 

2Cj 

AP*/'T* 

0 

0 

1.38 

1.40 

0.350 

- 

0.372 

0.222 

0.2 

1.37 

1.39 

0.352 

0.178 

0 . 504 

0.253 

0.333 

0.3 

1.35 

1.37 

0.354 

0.176 

0.500 

0.189 

0.444 

0.4 

1.32 

1.39 

0 , 360 

0.170 

0.489 

0.121 

c ,  660 

- 

- 

- 

- 

0 
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FIGURE  20 


DESIGN  AND  OPERATING  CHARACTERISTIC  CURVES 

The  graphical  results  for  the  sample  design 
calculations . 


FIGURE  21 


MAXIMUM  EJECTOR  EFFICIENCY  AS  A  FUNCTION  OF  GEOMETRY 


APPENDIX  A 

VELOCITY  AND  TEMPERATURE  TRAVERSE  RESULTS 


Velocity  aiid  Temperature  Profiles  --  Some  typical  ejector 
velocity  and  temperature  profiles  are  presented  in  Figs.  22  and 
23.  Figure  22  illustrates  the  profile  development  at  stations 
along  the  mixing  tube;  the  data  points  were  obtained  from  av¬ 
erages  on  either  side  of  the  tube  centerline.  The  x/D  =  0 
profiles  approximate  quite  closely  the  Idealised  entering  pro¬ 
files.  The  temperature  profile  at  x/D  =0  is  less  sharply 
divided  into  two  regions  due  to  heat  transfer  from  the  primary 
stream  through  the  copper  nozzle  to  tne  secondary  stream.  Al¬ 
though  this  effect  is  more  noticeable  at  low  flow  rate  ratios, 
the  influence  on  Km  is  always  small,  the  measured  being 

within  5  per  cent  of  the  values  reported.  The  presence  of  In¬ 
flection  points  indicates  that  the  profiles  are  far  from  es¬ 
tablished  at  x/D  =  h  ,  a  result  verified  by  the  high  K  at 
■this  station  (Fig.  16). 

Figure  23  describes  the  good  flow  symmetry  found  in  the 
ejector  at  an  x/D  of  2  by  superimposing  the  data  points 
taken  on  either  side  of  the  tube  centerline.  The  asymmetry, 
caused  by  nonuniform  entrainment  and  mixing,  is  seen  to  be 
negligible  for  the  temperature  distribution;  the  disagreement 
between  temperatures  on  either  side  of  the  centerline  is  about 
one  per  cent  of  full  scale.  The  velocity  distribution  is  more 
asymmetric;  the  difference  is  about  five  per  cent  of  full  scale. 
These  profiles  represent  the  maximum  asymmetry,  occurring  at 
low  w*T*0*^  and  an  x/D  of  2.  The  profiles  also  Illustrate 
the  flow  separation  at  the  wall  occurring  for  large  L./D  at 

x/D's  of  1,  2,  and  3  for  low  W*T*0*^;  U/b'  is  negative  for 
.  5s 

0.9  <  r/r0  <  l.C  .  The  separation  region  is  stable,  and  occu¬ 
pies  up  to  20  per  cent  of  the  mJ xing  tube  cross-sectional  area. 

As  stated  previously,  separation  is  not  present  for 
x/D=i  L/D  =  1,  2,  or  3  . 

Diffusion  of  Momentum  and  Temperature  --An  examination 
of  the  temperature  and  velocity  profiles  reveals  that  temperature 
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spreads  faster  than  momentum  (Pigs.  22  and  23).  This  conclu¬ 
sion  is  confirmed  by  the  work  of  Landis  and  Shapiro  (9);  how¬ 
ever,  the  statement  that  the  profiles  are  similar  is  not  con¬ 
firmed.  In  Fig.  23,  the  temperature  is  nearly  unifor;  across 
the  tube  while  the  velocity  variation  is  great  ar.d  separation 
occurs  at  the  wall.  This  disagreement  results  from  a  major 
difference  in  system  geometry.  The  results  of  (9)  were  taken 
on  an  ejector  with  an  area  ratio  of  63,  while  the  present 
ejector  has  an  area  ratio  of  1.835.  The  results  of  (9)  are 
to  be  interpreted  as  approximating  the  mixing  characteristics 
of  a  jet  discharging  into  an  infinite  flowing  medium,  and  not 
those  of  low  area  ratio  ejectors. 

The  centerline  decay  results  (Fig.  24)  are  compatible 
with  the  experimental  results  of  (9),  curves  b)  and  c)  , 
and  with  the  established  correlation  for  free  jet  discharge 
into  a  medium  at  rest,  curve  a)  .  The  trend  is  one  of  increas¬ 
ing  decay  of  centerline  velocity  wich  decreasing  w*T*0''  . 

The  established  empirical  correlation  for  centerline  velocity 
decay  for  jets  discharging  into  an  Infinite  flowing  medium  is 
given  in  (10)  as: 

ue  -  V  L 

-± -  =  -  ,  for  x/D  >  L 

VP  -  Vs  */D„ 


where  L 
by 

and 


the  length  of  the  potential  core  region,  is  given 

L  =  4  +  12A 

A  =  V  /V 
s  p 

x/D  =  flow  length/primary  nozzle  diameter 

*  j 


Attempts  to  use  this  correlation  to  predict  the  present  ejec¬ 
tor  behavior  are  also  shown  on  Fig.  24.  The  agreement  between 
the  data  points  for  x/Dn  -  6.6  and  the  prediction  indicates 
that  the  correlation,  for  T*  =  1.00  ,  may  hold  for  low  area 
ratio  ejectors.  However,  os  data  points  with  different  A  , 
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but  equal  W*T*  ,  have  the  same  centerline  decay,  A  Is 
r.ot  the  only  parameter  governing,  the  centerline  velocity  de¬ 
cay.  In  terms  of  the  present  nomenclature,  A  can  be  ex¬ 
pressed  as: 


or 


w 


PsAs 


A  p  W*T* 

_£kP_  ^  _ 

w  A* 

s 


,  v0.5  W*T*°*5 

A/  I*  =  - r* - 


In  view  of  the  present  experimental  results,  it  seems  prob¬ 
able  that  A/T*^*^  ,  rather  than  A  ,  should  be  used  in  the 
correlation.  However,  more  data  is  needed  to  definitely  con¬ 
firm  this  choice. 

The  centerline  temperature  decay  results.  Fig.  25,  do 

not  follow  the  trends  of  (9)  ss  closely  as  the  velocity  decay 

results  do.  The  behavior  is  similar  in  that  the  decay  in- 

0  5 

creases  as  the  value  of  A  ,  and  W*T*  decreases,  but  the 
lack  of  enough  data  precludes  a  meaningful  comparison.  Th= 
difference  in  geometry  influences  the  centerline  decay  also. 

The  existing  centerline  decay  correlations,  such  as  from 
(9),  are  derived  mainly  from  high  area  ratio  ejectors  and  are 
not  suitable  for  low  area  ratio  ejectors  as  these  correlations 
require  that  the  centerline  values  approach  the  secondary 
stream  values  as  x/C^  becomes  very  large.  For  high  area 
ratio  ejectors,  this  is  a  fairly  good  approximation  as  can 
be  demonstrated  by  the  following  considerations.  For  all  ejec¬ 
tors,  the  flow  at  large  x/D^  approaches  fully  established 
turbulent  pipe  flow  and  the  centerline  velocity  approaches 
within  20  per  cent  the  bulk  average  velocity  (U„  >  V  always). 

Vm  is  essentially  equal  to  Vg  if  the  primary  area  is  small 
compared  to  the  secondary  area  (large  A*  ).  Also  for  fully 
established,  low  velocity,  aaiabatic  flow,  the  centerline  tem¬ 
perature  is  equal  to  the  bulk  average  temperature,  and  Tm  is 
essentially  equal  to  Tg  if  A*  is  large  so  that  the  primary 
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flew  rate  is  small  compared  to  the  secondary  flow  rate.  For 
low  area  ratios,  and  therefore  low  flow  rate  ratios,  this  is 
not  true;  the  bulk  average  quantities  are  considerably  differ¬ 
ent  from  the  secondary  quantities,  and  the  temperature  and 
velocity  decay  parameters  will  not  tend  to  zero  as  x/Dn  in¬ 
creases.  This  can  easily  be  seen  by  formulating  the  bulk  av¬ 
erage  velocity  and  temperature,  normalized  by  V  and  T 

Jr  J. 

respectively,  in  terms  of  the  present  nomenclature; 


V 

m 

V 

P 


(w  +  w  )  p  A  (1  +  W*) 
2  p  sy  _  *p  p  _ _ _ 

(A  +  A  )p  v/  (l  +  A*) 
'  p  s  '  ■  m  p  '  ' 


^m  _  wp  #  +  ws 

Tp  <wp  +  ws>  Tp  <wp  +  ws> 


(1  +  W*T* ) 
(1  +  W*) 


The  limiting  values  of  the  nondimensional  parameters  used  to 
correlate  the  centerline  decays  can  now  be  formed.  The  fully 
established  value  of  Ug  is  assumed  to  be  1.2Vm  ;  the  fully 
established  value  of  T„  is  T 


-  V 


V  -  V 
P  s 


1  (1  +  W*T* ) 

-  1  +  0.2A*  - 

1  +  A*  (A*  -  W*T* ) 


T„  -  T 


T  -  T 
P  s 


1  +  W* 


For  the  large  area  ratio  ejector,  as  in  (9),  these  limits 
are  relatively  close  to  ze~o;  for  A*  =  65  and  W*T*^ * ^  =  56  , 


U„  -  V. 


V  -  V 
P  s 


=  0.28 


T.,  -  T 

_L - d  =  0.027 

T  -  T 
P  ~s 


However,  for  low  area  ratio  electors,  these  limits  are  con¬ 
siderably  different  from  zero;  for  A*  =  1.855  and 
W*T*0,5  =0.64  , 
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and 


0.53 


=  o.6l 


The  conclusions  and  recommendations  that  can  be  drawn 
from  this  discussion  are: 


1.  Further  investigations  of  the  temperature  and  mo¬ 
mentum  profiles  are  needed  to  establish  the  mixing  profiles 
of  low  area  ratio  ejectors  and  to  establish  the  correspond¬ 
ence  between  low  and  high  area  ratio  ejectors. 

2.  More  data  is  needed  to  ascertain  if  the  empirical 

correlation  presen:ed  in  (10)  for  the  centerline  velocity 

decay  is  valid  for  low  area  ratio  ejectors  for  even  a  limited 

0  *5 

region,  and  whether  A/T*  is  a  better  correlating  parame¬ 
ter  than  A  for  centerline  decay.  It  is  also  necessary  to 
investigate  the  importance  of  the  bulk  average  properties, 
and  to  determine  whether  these  properties  can  profitably  be 
used  in  establishing  a  correlation. 


To  this  end  it  is  recommended  that  data  be  taken  of  the 
centerline  velocity  and  temperatures  over  a  range  of  A*  and 
T*  magnitudes.  These  results  should  provide  a  test  for  the 
validity  of  the  proposed  A/T*  correlating  parameter.  It 
also  seems  probable  that  the  velocity  decay  should  be  expressed 
as 


(UE  -  vm> 


<7p  -  V 


rather  than 


-  vs> 

-  V 


and  similar  ratios  for  temperature.  Then  both  the  low  and 
high  A*  geometries  would  have  common  asymptotic  magnitudes 
for  large  x/Dn  . 
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FIGURE  22 


VELOCITY  AND  TEMPERATURE  PROFILES 

Illustrates  development  of  the  profiles  at 
stations  along  the  mixing  tube. 


FIGURE  23 

VELOCITY  AND  TEMPERATURE  PROFILES  AT  A  LOW  x/D  AND  W*T*°‘5 
Illustrates  the  asymmetry  ana  separation  present. 


FIGURE  24 

centerline  velocity  decay 


FIGURE  25 

CENTERLINE  TEMPERATURE  DECAY 
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AfPiiNDlX  B 

ANALYSIS  OF  EXPERIMENTAL  UNCERTAINTIES 

The  estimated  experimental  uncertainties  of  the  main 
test  results  are  as  follows; 


Quantity 

Uncertainty 

W* 

i  1.2$ 

T* 

±  0.1$ 

W*T*°.5  and  w*t*0*4 

±  1.2$ 

AP*/T* 

±  3.1$ 

Km 

(from  a?*/t*) 

±1  $ 

Km 

(x/D  >  4  ,  traverse ) 

±2  $ 

Km 

(x/D  <  4  ,  traverse) 

+1  to  +5$ 

The  Xm;s  evaluated  from  the  velocity  and  temperature 
traverses  have  higher  uncertainties  compared  to  the  evalua¬ 
tions  from  overall  ejector  performance  as  the  highly  turbu¬ 
lent  and  rotational  flow  in  the  mixing  regions 'produce  large 
uncertainties  in  the  profiles  (a  discussion  of  these  effects 
follows  in  Appendix  C),  The  uncertainty  due  to  these  fluc¬ 
tuations  is  determined  by  the  following  observations. 

The  pitot  tube  readings  are  too  high  due  to  turbulent 
velocity  fluctuations;  the  velocity  profiles  computed  from 
these  readings  will  be  too  high  also.  As  the  velocity  gra¬ 
dient  is  greatest  between  x/D  =  0  and  x/P  -  4  ,  the  fluc¬ 
tuations  will  be  highest  here.  This  is  evidenced  in  that 
the  discrepancy  between  metered  and  Integrated  flow  rates  is 
a  maximum  at  x/D  =  2  ,  and  negligible  at  x/D  =  6  and  8. 

A  crude  analysis  of  the  effect  of  these  fluctuations  was 
undertaken  on  Run  51  (W*  =  0.75,  x/D  =  2)>  Run  64  (W*  -  0.10, 
x/D  =  2),  and  Run  49  (W*  =  0.75,  x/D  -  4).  Corrisin  (4)  has 
measured  values  of  turbulence  in  free  jet  mixing;  from  this 
report,  fairly  typical  values  of  u'/Ur  =  0.50  and 
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v'/U  =  0.25  were  assumed  for  the  mixing  regions.  An  ex¬ 
amination  of  the  velocity  profiles  reveals  a  definite  mixing 
region  0.2  <  (r/r  )2  <  0.7  for  Run  51  (Fig.  26), 


0-1  <  (r/r0) 


<  (r/rQ)  <  0.7  for  Run 
—  0.9  U, 


for  Run  64,  and  u 

49.  The  velocity  correction,  Uactual - *  “indicated 

(see  Appendix  C  for  derivation),  was  applied  to  the  profiles 
in  these  regions.  The  values  of  flow  rate  and  were 

found  to  be  corrected  by  the  following: 


Run 

51 

64 


x/D 

W  /  w 

m,corr/  m,uncorr. 

K  /K 

m,corr/  m,uncorr. 

Increase 
in  K 

m 

A 

0 . 945 

1.034 

3.40 

2 

0.918 

1.028 

2.80 

h 

*T 

0.998 

1.003 

0.30 

As  the  uncertainty  in  determining  the  Km's  is  ±  20, 
the  total  uncertainty  in  Km  due  to  both  Instrument  uncer¬ 
tainty  and  turbulent  velocity  fluctuation  is  taken  to  be: 


x/D  >  4  ±  2% 

x/D  <4  +1  to  +50 
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FIGURE  26 


MOMENTUM  AND  FLOW  RATE  PROFILES 

Illustrates  effect  of  turbulence  on  thu 
mixing  tu1. e  flow  rate  and  on  . 


APPENDIX  C 

ANALYSIS  OF  PITOT  TUBE  RESULTS 


One  of  the  primary  objectives  of  this  work  was  to  es¬ 
tablish  confidence  in  the  pitot  tube  traverse  results  by  a 
comparison  of  the  metered  flow  rate  and  that  obtained  by  a 
velocity  traverse.  It  was  observed  that  tne  traverse  re¬ 
sults  were  consistently  high.,  especially  at  low  x/D's  - 
from  2%  at  a  W*  of  1.0  to  10$  at  a  W*  of  0.1.  This  sug¬ 
gested  a  re-examination  of  Pitot  tube  results. 

Goldstein  (6)  establishes  that  the  equation  of  motion 
in  the  longitudinal  direction  for  turbulent  flow  is 

P  1  tt2  ,  1  —  2  Pt 
p+2U+2q  -  p 

* 

where  P  is  the  free  stream  pressure;  U,  the  flow  mean  ve¬ 
locity;  and  q"1  ,  the  resultant  turbulent  velocity  fluctua¬ 
tion, 

—  —  2  —  P  —  2 

q'  =  U'  +  v*  +  w' 

The  equation  of  motion  in  the  transverse  direction  is 

P  1  „  Pc 

p  +  2  v  '  =  p~ 

where  v'  Is  the  transverse  velocity  fluctuation.  The  pitot 
static  tube  equation  then  becomes 

(rt  -  p„)  l  ,  i  _  -  _  , 

— - —  =  -  u2  +  -  (q'2  -  v'2) 

p  2  2 

or  U  ,  the  mean  velocity,  must  be  calculated  from 


A  /  (Pt  Ps)  ^  2  —  2\ 

=  Vs  ~  (q'2  -  V2) 


In  most  turbulent  flows,  the  turbulent  components  are 
negligible.  However  Corrlsin  (4)  has  studied  the  turbulent 
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fluctuations  in  free-jet  mixing  and  found  values  of  the  ratio 
of  the  turbulent  velocity  fluctuation  to  local  velocity  in  the 
neighborhood  of  U'/U  =  0.5C  and  "v'/dr  =  0.25  for  the 
highly  turbulent  regions.  These  values  of  the  fluctuating 
components  would  result  in  a  pitot-static  reading 


1.19  U2/2 


The  error  .in  velocity,  if  the  turbulent  fluctuation  is  neglected 
in  the  calculation,  is  9$.  The  effect  on  the  traverse  results 
by  neglecting  the  turbulent  fluctuations  has  also  been  observed 
by  Nielsen  (7)  in  comparing  the  flow  metered  by  an  orifice  to 
that  obtained  by  the  Integration  of  the  velocity  profile,  the 
traverse  being  made  a  small  distance  behind  the  orifice.  Here 
flow  rate  discrepancies  from  >  to  16$  were  observed. 

For  highly  turbulent  regions  of  flow,  some  knowledge  about 
the  turbulent  fluctuations  is  required  if  the  velocities  and 
flow  rates  are  to  be  determined  from  Pitot  traverses.  If  the 
data  are  not  corrected  for  this  effect,  the  calculated  veloc¬ 
ity  will  be  higher  than  the  actual,  and  the’  flow  rate,  propor¬ 
tional  to  the  velocity,  will  likewise  be  high.  This  effect 
also  Indicates  more  complete  mixing  than  actually  occurs,  with 
the  result  that  the  calculated  K  's  are  lower  than  the  actual. 
These  conclusions  are  supported  in  Appendix  B,  where  a  veloc¬ 
ity  profile  was  adjusted  for  this  effect. 

Another  possibility  for  discrepancy  exists  if  there  is 
an  appreciable  transverse  velocity  component.  In  the  mixing 
tube,  this  could  be  flow  that  spirals  down  the  tube.  A  dis¬ 
crepancy  is  caused  in  that  the  pitot  tube  measures  the  total 
velocity  head,  even  if  misaligned  as  much  as  20°  with -respect 
to  the  mean  flow.  The  resultant  velocity,  of  greater  magni¬ 
tude  than  the  longitudinal  or  mean  velocity,  will  be  interpre¬ 
ted  as  the  longitudinal  velocity.  The  calculated  velocity 
will  then  be  larger  than  the  actual  longitudinal  velocity. 

The  flow  ang’e  can  be  determined  with  a  yaw  detector.  The 
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results  Indicated  an  angle  of  about  2.5°,  or  0.1#  on  veloc¬ 
ity.  The  effect  was  negligible  in  this  investigation. 

A  third  type  of  error  is  introduced  into  the  pitot  cube 
results  in  that  the  flow  is  highly  rotational.  This  causes  a 
shifting,  or  curving, of  the  streamlines  in  the  vicinity  of  the 
pitot  tube  nose  such  that  the  pitot  tube  "sees"  a  velocity 
higher  than  the  velocity  that  would  be  there  in  the  absence 
of  the  pitot  tube.  An  empirical  relation  for  the  amount  of 
shift  ,1s  given  by 


D/Dto  =  0.131  +  0.C82  Dt0/Dtl 


D.n  dU 

0.1  <  -12. - <  1.2  ,  (8) 

U  dr 


where  DtQ  is  the  pitot  tube  outer  diameter,. 


0.058";  Dtj  , 


the  inner  diameter,  0.029"  ;  (dU/dr)  ,  the  velocity  gradi¬ 


ent;  and  6  ,  the  shift  toward  the  low  velocity  region, 

6  —  0.01"  ,  o/rQ  —  0.007  •  As  the  flow  at  x/D  <  2  was  of 
this  type,  this  correction  was  effected  on  the  raw  data." 
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